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ABSTRACT 

We present the complete data set, model and line identification of a survey of the emission from 
the C-rich protoplanetary nebula CRL 618 performed with the IRAM-30m telescope in the following 
frequency ranges: 80.25-115.75 GHz, 131.25-179.25 GHz, and 204.25-275.250 GHz. A selection of lines 
from different species has been used in previous works to derive the structure of the source, its physical 
conditions and the chemical abundances in the different gas regions. In this work, we have used this 
information to run a global simulation of the spectrum in order to check the consistency of the model 
and to ease the task of line identification. The total number of lines that have a correspondence in 
both data and model is ~3100, although quite often in this object many lines blend into complex 
features so that the model, that takes into account line blending, is a key tool at this stage of the 
analysis. Of all the lines that we have been able to label, ^55% of them belong to the different forms 
of HC3N, and ~18% to those of HC5N. The density of remaining unidentified features above the 3cr 
limit is only one per ^2.1 GHz (74 features), which is unprecedented in the analysis of this type of 
large millimeter- wave line surveys. 

Subject headings: line identification - surveys - stars: post-AGB-stars: carbon-rich - stars: circumstcl- 
lar matter - stars: individual: CRL 618 - ISM: molecules - radio lines: stars 
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1. INTRODUCTION 

In order to study the evolutionary stages from Asymp- 
totic Giant Branch (AGB) to Planetary Nebulae (PN), 
multiwavelength and multi-object data are necessary to 
reveal the key changes in physical conditions and chem- 
ical composition as completely as possible. Even so, the 
task is not trivial due to the many distinct environments 
that can be present in these C-rich or O-rich objects in- 
cluding: very fast winds and presence of shocks (Cer- 
nicharo et al. 1989; Davis et al. 2005; Bujarrabal et 
al., 2002), binarity (Sanchez-Contreras et al. 2004, De 
Marco et al. 2004), magnetic fields (Frank & Blackman 
2004; Huggins k. Manley 2005), late thermal pulses (As- 
plund et al. 1999, Pavlenko et al. 2004), etc... Molecular 
spectroscopy at millimeter, submillimeter and FIR wave- 
lengths is a useful tool to study the early stages of this 
evolution because large carbon chain and metal contain- 
ing molecules form in the atmospheres of late-type stars. 
This rich formation chemistry exists in C-rich objects due 
to the relatively high densities and the presence of dust 
and shocks created by significant mass loss (Cernicharo 
2004). In addition, the potentially photodissociating ra- 
diation from the central star, when it changes its spectral 
type, is blocked in a compact inner region. As the evolu- 
tion goes on, the shielding of high energy photons from 
the central star, that has evolved toward types A, B and 
O is less effective so that the molecular content in the 
envelope decreases greatly. 
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The ideal targets to study the formation chemistry 
of large carbon-chain molecules are late type stars with 
significant mass loss (e.g., IRC+10216), protoplanetary 
nebulae (e.g., CRL 2688 and CRL 618) and very young 
planetary nebulae (e.g., NGC 7027). Herpin et al. 
(2001), performed a study using several CO lines across 
the mm and submm ranges in the latter four objects, 
taken as prototypes for the different steps of the AGB to 
PN evolution, aimed at determining the physical condi- 
tions prevailing for the molecular gas. Sanchez-Contreras 
et al. (2002, 2004) have undertaken a multiwavelength 
study of CRL 618. Cernicharo et al. (1996, 2001a,b) have 
explored the ISO data of CRL 618 and IRC+10216 with 
interesting findings such as the first evidence for the first 
aromatic molecule seen outside the solar system (CgHg 
in CRL 618). Remijan et al. (2005) have performed a 
search for large molecules toward CRL 618 although no 
detection of biologically important molecules could be 
reported. Finally, Cernicharo et al. (2000) published 
a A=2 mm survey of IRC+10216. Therefore, our mo- 
tivation for the complete Institute de Radioastronomie 
Millimetrique (IRAM) 30 meter telescope survey of CRL 
618, fully presented here, has been to gather the most 
complete molecular information on one particular stage 
(protoplanetary nebula, PPNe) of low-mass stellar evolu- 
tion, of which CRL 618 is the best C-rich example. The 
data on their own have allowed a quite precise determi- 
nation of physical conditions and molecular abundances 
in this object (Pardo et al. 2004, 2005; Pardo & Cer- 
nicharo 2006; hereafter P04, P05 and PC06). These pre- 
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vious works have allowed us to build a complete model 
to be compared with the whole data set, thus simplifying 
the task of isolating U-lines during the line identification 
process. 

2. OBSERVATIONS SUMMARY 

The CRL 618 IRAM-30m line survey was completed in 
1994-2002. In order to keep the calibration of the spec- 
tra as consistent as possible, the observing procedure was 
kept unchanged through the years. Nevertheless, techni- 
cal works at the telescope and varying atmospheric condi- 
tions are responsible for a relative calibration uncertainty 
of the order of 10%. This uncertainty in the calibration 
was based on observing the same lines at different epochs. 
The analysis of more than 15 lines in different vibrational 
states of HC3N has allowed to notice a few scans with ob- 
vious calibration problems (less than 10 over more than 
350 spectra, or <3%) most probably due to an erroneous 
calibration scan. Those few spectra were re-scaled based 
on the overall results of the rotational HC3N ladders. 

In order to avoid signal and image sideband confusion, 
image sideband rejection larger than 12 dB was used. 
Features that could be attributed to sideband contami- 
nation are only a few, generally related to a very intense 
line (CO, HCN, HNC, HCO+,...) in the image side band. 
The pointing and focus were always checked on the tar- 
get source and the peak level of the gaussians after the 
last correction of each pointing session has been used to 
study the spectral behavior of the millimeterwave con- 
tinuum in CRL 618. The pointing itself was kept within 
2" accuracy. The observations were carried out using the 
wobbler-switching mode with offsets of 60" and frequen- 
cies of 1 Hz, in order to obtain very flat baselines. The 
backends were two 512x1 MHz filter banks connected to 
the receivers operating below 200 GHz, and an autocor- 
relator with channel widths of 1.25 MHz connected to 
the receivers operating in the 1.3 mm window. System 
temperatures were typically 100-400 K at 3 mm, 200-600 
K at 2 mm, and 300-800 K at 1.3 mm. 

Complementary observations above 280 GHz have 
been carried out with the 10.4 m dish of the Caltech Sub- 
millimeter Observatory (Mauna Kea, Hawaii) and some 
of them have been used in the analysis presented in P04, 
P05 and PC06. The complete results above 280 GHz will 
be published elsewhere. 

3. RESULTS 

In order to illustrate the overall results of the survey we 
include here two tables (0]and[5]) describing the molecules 
detected and the J up or energy range in which the lines 
of a particular species are seen. Other species included in 
the model, although marginally or not detected, appear 
also in the same table. We also provide four figures show- 
ing the overall data and model in four spectral ranges, 
corresponding to those of the available receivers at the 
IRAM-30m telescope: 3 mm, 2 mm, 1.3 mm (lower end) 
and 1.3 mm (upper end). In order to have an idea in 
the printed version of the correspondence between data 
and model, we present nesting zooms in frequency in all 
figures, so that the most detailed panel shows 0.5 GHz 
wide spectra (corresponding to the frequency coverage 
in each observational setting). The whole survey at this 
level of detail needs about 300 pages of figures and ta- 
bles and therefore this is available only as on-line mate- 



rial. Finally, Table [6] provides the position and intensity 
of the remaining unidentified features after the analysis 
performed. 

3.1. Continuum emission 

The continuum in T^ scale was derived from the point- 
ing scans, as explained in section[2l after discarding a few 
obvious bad scans. Tj\/s and fluxes have been derived 
from a fit of those data to a horizontal line, scaled by the 
corresponding beam and forward efficiency factors of the 
30m telescope (B e // ranges from 0.79 to 0.42 and F e ff 
from 0.95 to 0.88 in the frequency range 80 to 275 GHz 
for this telescope). This analysis is not too much affected 
by the presence of lines because the total flux from the 
lines is less than 3-5 % relative to the continuum flux 
in most frequency settings. As a comparison, the total 
flux is dominated in fact by the lines in millimeter wave 
observations toward high mass star forming regions such 
as Orion (> 50% of the total flux towards the IRc2 po- 
sition at these wavelengths, Tercero et al. in prep). A 
few frequency settings with very strong lines such as the 
lowest rotational transitions of CO and HCN have been 
ignored for the continuum flux analysis in CRL 618. The 
results were presented and discussed in sections 4.1 and 
5.1, Figl, and Table 2 of P04. In terms of T M s, the 
continuum follows a straight line from 0.52 K at 80 GHz 
to 1.52 K at 275 GHz. The spectral behavior seems to be 
basically the same in the frequency range 280-360 GHz 
according to our results obtained at the Caltech Submil- 
limeter Observatory (P04, P05). 

3.2. Overall line spectra 

The overall results of the line survey presented in fig- 
urcs[TJ[2l[3l and[4]show that the most intense lines (above 
0.2 K in T^) belong to a very limited number of species 
(H recombination lines, C 3 H 2 , HCN, HCO+, CN, CO, 
CS, H 2 CO, CH 3 CN, CCH and HC 3 N). They also show 
how, as frequency increases, the absorption part of the 
line profile becomes less important with respect to the 
emission part in species showing P-Cygni profiles, indi- 
cating that the spectral behavior of the continuum flux 
plays a fundamental role in the evolution of the absorp- 
tion/emission ratio. The absolute energy levels of the 
transitions play also a role, but less important. The 
clustering of lines belonging to HC3N and its isotopic 
and vibrationally excited species every ~9 GHz is also 
nicely delineated. The detailed zoom (500 MHz) at the 
bottom of each figure illustrates the line identification 
process. Individual figures at this level of detail, includ- 
ing a comparison with the model, have been created for 
the whole data set and are available only as on-line mate- 
rial. Some spectra are so crowded with lines that zooms 
showing only 250 MHz of data have been necessary. A 
set of tables with the format of Table ??, only available 
also as on-line material, provide the rest frequencies of all 
identified features and the complete quantum numbers of 
the transition (usually omitted in the figures due to lack 
of space). The labeling is quite straightforward for all 
molecules with the exception of HC3N that, we should re- 
mind, is responsible for about 55% of the features seen in 
the survey. For this molecule the lines are labeled by the 
rotational quantum numbers (J M p-Ji ou) ), the vibrational 
quantum numbers of the four (out of seven) lowest vi- 
brational modes {yi,V5,v§,v>j) and the ^-doubling param- 
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eters for the bending modes: [£5 £g £7 ]. The complete 
set of frequencies for this molecule has been obtained 
from Fayt et al. (2004), and the agreement is excellent. 
For HC5N we use a similar notation with quantum num- 
bers of the three (out of eleven) lowest vibrational modes 
(vg,vio,vn) and the necessary ^-doubling parameters. 

3.3. Line identification 

The line identification process has been based on J. 
Cernicharo's own molecular rotational line catalog. The 
catalog contains now ~1300 molecules, radicals and 
atoms, and the number of lines in the frequency range 
80-280 GHz exceeds 600000. 

Although part of the line identification could be done 
by means of an automated procedure (mainly designed 
to recognize the lines from the cyanopolyynes HC3N and 
HC5N), the spectra are sometimes so crowded with lines 
that an important part of the identification and label- 
ing work has to be done manually, specially taking into 
account that the greatest interest is usually in the weak 
features. 

The molecular species with identified lines in the mil- 
limeter spectrum of CRL 618 are listed in Tables 0] and 
[3 Very few are O-bearing molecules (only SiO, H2CO 
and HCO + in addition to CO), and there is also CS, 
and recombination lines of Hydrogen and Helium. The 
most prominent features in the spectrum are due to the 
cyanopolyynes family, mainly HC3N, for which we have 
identified lines from more than 30 isotopic, isomeric or 
vibrationally excited forms. The key for deducing the 
basic physical parameters of CRL 618 has been to count 
on several hundred of lines of HC3N with a wide range of 
excitation conditions to probe the different gas regions 
of this source (see P04 and P05). 

We know nevertheless that this object has to have some 
relatively large carbon chains, based on the discoveries 
published by Cernicharo et al. (2001a, b). Compared to 
other cases such as Orion or IRC-l-10216, the number 
of U- lines in the CRL 618 survey is very small: Only 
74 or one per ~2.1 GHz as an average in the surveyed 
range. Although several species are candidates to ex- 
plain a few of these U lines, their complex rotational 
spectrum should display other features arising above the 
noise level that are not seen, and therefore the assign- 
ment could not be established. Table [6] gives the fre- 
quencies and peak emission and absorption (if present) 
values for the unidentified features after subtracting the 
continuum level. The search for some particular species 
is described in section l3~5l 

Detailed tables and figures showing the ^3100 features 
that could be identified with the help of the model de- 
scribed in section 13.41 are available as on-line material 
only. 

3.4. Model 

A model aimed at reproducing the whole millimeter 
wave spectrum of CRL 618 as accurately as possible with 
the minimum number of parameters has been developed. 
To achieve its completion, four basic steps have been 
followed: 

1. Model for the continuum flux considering a fixed 
size of the region from which it emerges, its bright- 
ness temperature at a given frequency, and a spec- 



tral index. The values used, 0.27" size, 3900 K at 
200 GHz and (y /^o) -1 ' 12 behavior, were found and 
discussed in P04. 

2. Model for the inner (1.5") slowly expanding molec- 
ular envelope (SEE) using lines from vibrationally 
excited states of HC3N (the best tool for this task). 
The average temperature of the region (250-275 K), 
HC3N column density in front of the continuum 
source (~3T0 17 cm -2 ), and details about the mor- 
phology and the velocity field, were also obtained 
in P04. 

3. Model for the high velocity wind (HVW), reaching 
^200 kms -1 , presumably quite collimated, notice- 
able in some abundant species such as CO, HCN 
and v=0 HC 3 N (see P05). 

4. Model for the colder and outer, with respect to the 
SEE, cold circumstellar shell ( CCS) that is respon- 
sible for most of the rotational emission from v-i 
and v=0 HC 3 N, and of v=0 HC 5 N. This gas com- 
ponent is also dominant in the emission of several 
other species. The average temperature (~30 K), 
position (3.0" to 4.5" diameter), velocity field, den- 
sity distribution, and abundance ratio of all species 
with respect to HC3N, have been obtained in P05 
and PC06. 

The model goes into some detail on the morphology of 
the object. For example, the geometry is not spherical. 
The SEE has an elongated shape, the HVWis collimated, 
and the CCS has a bipolar density distribution. The ve- 
locity field has both radial and azimuthal components. 
The azimuthal symmetry is broken because the line of 
sight is inclined with respect to the symmetry axis. See 
Fig. 4 of P04 and Fig. 3 of P05 for more details. As a re- 
sult, as discussed in section 3.2 of PC06, quadratures in 4 
different variables are necessary for the calculations (fre- 
quency, impact parameter, path along the line of sight, 
and angle in the plane perpendicular to the line of sight). 
All species listed in Table 0] have been included in the 
general calculations (except stated otherwise). 

In P04, P05 and PC06, the calculations were performed 
for a selection of lines of each species that met several 
requirements such as avoiding line blendings and repre- 
senting well the energy range probed by each species. 
This allowed us to draw a quite precise picture of the 
structure, physical conditions and chemical abundances 
in the source with reasonable computation times. Now, 
the model has been run for all surveyed frequencies in 
order to check its consistency and to be provided as a 
final product of this study together with the whole and 
labeled data set. 

3.5. Search for Trace Molecules 

As the complexity of the molecules increases, their 
abundances are usually lower and their partition func- 
tions larger. In order to identify a new molecule in the 
data, it is necessary to run a model for a large set of tran- 
sitions. Our search for new species has been based on 
possible products from simple reactions using molecules 
and radicals already identified in this survey and in IR 
observations. 
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3.5.1. Ce-ff radical 

The abundance of C2H and CN in CRL 618 allows an 
efficient growth of HC2 n +iN and C271H2 molecules (see 
the detection of the polyyne chains C2H2, C4H2, and 
CgH2 in Cernicharo et al. 2001; and the polymerization 
of HCN in Cernicharo 2004 and P05). As a consequence, 
the linear radicals C2„H, with only one free bond, are 
much more abundant in this object than C2 ra +iH, with 
three free bonds. Therefore, we expected to detect the 
CqH linear radical. However, the large rotational parti- 
tion function of this species combined with a decrease in 
abundance with respect to C4H and also the the higher 
energy levels of CgH that are sounded above 80 GHz (well 
above the Boltzmann peak corresponding to a tempera- 
ture of ~30 K) with respect to the lighter radicals, make 
in fact its detectability marginal at the sensitivity level 
of the survey. 

3.5.2. C S N radical 

Although the radical C3N is more abundant than the 
combination of C-C3H and I-C3H, C5N is not detected 
in the survey, contrary to C5H. The latter is the less 
abundant main species detected in the survey. Assum- 
ing a similar abundance for C5N its non detection could 
be explained by its lower dipole moment (3.385 Debyes 
instead of 4.881 Debyes for C5H). In any case, above 80 
GHz the rotational lines of both molecules are past the 
Boltzmann peak corresponding to the temperature of the 
cold circumstellar shell of CRL 618. 

3.5.3. H2C4 and H2C3 radicals 

Double C bounds seem to be rare in CRL 618 com- 
pared with carbon chains alternating single and triple 
bounds (as in the cyanopolyynes, polyynes and related 
radicals). This is further confirmed by our search of the 
H2C4 and H2C3 radicals (in which all carbons are dou- 
bly bounded) that has resulted on non detection for H2C4 
and a marginal one for H2C3. Although some features in 
the spectra are located at H2C4 and H2C3 frequencies, 
a calculation of the whole spectrum to reproduce those 
features gives inconsistencies at many other frequencies 
(non-detections where a line should be expected). 

3.5.4. CH 2 CHCN 

This is so far the only positive new finding. This 
molecule has /i a =3.815 Debyes and /if,=0.894 Debyes 
transitions that result on many weak features in our sur- 
vey (see Fig. [21 bottom panel). We have found an abun- 
dance ~15 times less than that of HC3N in the CCS. This 
value is intermediate between those found for CH3CN 
and CH3CCH, although there is evidence that part of the 
signal detected from the last two molecules comes from 
the SEE (see PC06). The relatively similar abundances 
found for the three species tell us about a similar chem- 
istry leading to their formation; see Cernicharo (2004) 
for details about chemical routes in CRL 618. 

3.5.5. C 6 H 5 CN and C 6 H 5 CCH 

Given the evidence in support of the detection of ben- 
zene (CgHg) and the large abundances of the radicals CN 
and CCH in CRL 618, we have checked for the presence 
of these two molecules but unfortunately with negative 



results. The large partition function of these species com- 
bined with a column density of benzene of 2xl0 15 cm~ 2 
at most in this object (Cernicharo et al. 2001) in fact 
imply that the intensity of the lines cannot be expected 
to arise above the sensitivity level of the survey. 

3.5.6. CH 2 CN 

This radical, that in principle could be expected to 
be detectable as it should be formed in similar routes 
to those of the abundant CH3CN in CRL 618 has been 
the subject of an specific search. Only some features 
around 161.95 GHz and 241.35 GHz (see table [6]) seem 
to be related to this radical, but the detection cannot be 
confirmed because there is no correspondence at other 
frequencies where a line should be expected. 

3.5.7. 13 C substitutions of C-C3H2 

Starting from the abundance of C-C3H2 and consider- 
ing a ratio 12 C/ 13 C of 40 in the CCS, found in P05 and 
PC06, it would result that H 13 CCCH at least could be 
detectable. Its non detection indicates that perhaps the 
value 40 should be understood as a lower limit. Further 
evidence for this is provided by H^CO, 13 CS, 13 CCH, 
and C 13 CH, that are not seen either. We should remind 
that the ratio 12 C/ 13 C is only about 15 in the inner SEE, 
possibly due to injection of 13 C-rich material generated 
in a late thermal pulse to the inner parts of the gas out- 
flows of CRL 618. The result presented in this section, 
therefore, further strengths our (PC06) finding of quite 
different 12 C/ 13 C abundance ratios in the SEE and CCS 
regions of CRL 618. 

3.6. Comparison with PPN chemical models 

The study of Woods et al. (2003) provides chem- 
ical abundances in a region within 8.9T0 15 cm from 
the central star, which roughly coincides with the 
SEE of our model. Their predictions provide small 
[HC 3 N]/[HC 5 N]/[HC 7 N] ratios, just as we observed, 
that support the efficiency of polimerization in creating 
polyyne and cyanoplyyne chains in this inner SEE (al- 
ready reported in P05 and in agreement with predictions 
in Cernicharo 2004). The predictions made in Woods et 
al. (2003) for radicals such as CN, CCH, C 4 H, C 3 N,... 
show relative abundances with respect to HC3N of less 
than 10 -2 (except — 10 _1 for C3N) which would make 
them very difficult to detect in the SEE. This is con- 
firmed by our data because the emission we detect from 
those radicals is dominated by the CCS component (al- 
most no signs of absorption) . The reason for this is that 
at the densities and temperatures of the SEE, reactions 
with H2 would be very efficient to transfer those radicals 
into polyyne or cyanopolyyne species (that dominate in 
this region) whereas in the much colder CCS (T~30 K) 
the reactions with H2 would barely operate. 

4. CONCLUSIONS 

This paper has presented the final products of the most 
complete spectroscopic study at millimeter wavelengths 
of a protoplanetary nebula carried out up to date. First, 
the whole data set with line identifications for about 3100 
features has been presented. The number of remaining 
unidentified features is unprecedentedly low, only one per 
-2.1 GHz (a total of 74 in 155 GHz). Second, the whole 
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TABLE 1 

Identified features in the 80.3-115.5 GHz range in the CRL 618 survey presented in this paper. 



i/ (MHz) 


Molecule 


Transition / estate 


1/0 (MHz) 


Molecule 


Transition / t)-state 


80305.928 
80388.107 
80420.646 
80446.196 
80618.000 


HC 5 N 

i-c 3 h 2 n 3/2 
1-C3H 2 n 3/2 
hc 5 n 
U080618 


J=30-29 (003) [-1] 
(J,F)=(7/2,4)-(5/2,3) a 
(J,F)=(7/2,4)-(5/2,3) b 
J=30-29 (003) [1] 


80380.000 
80388.442 
80422.052 
80559.000 
80723.188 


HC5N 

1-C3H 2 n 3/2 
1-C3H 2 n 3/2 

U080559 
C 3 H 2 para 


J=30-29 (003) [-3] 
(J,F)=(7/2,3)-(5/2,2) a 
(J,F)=(7/2,3)-(5/2,2) b 

J/f+,/f-=42,2-4l,3 



Note. — Table Al is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 



TABLE 2 

Identified features in the 131.2-179.1 GHz range in the CRL 618 survey presented in this paper. 



i/ (MHz) 


Molecule 


Transition / u-state 


i/ (MHz) 


Molecule 


Transition / u-state 


131256.830 
131277.716 
131283.629 
131376.100 
132246.366 


MgNC 

c 5 h 2 n 1/2 

HC5N 
HC5N 
H 13 CCCN 


J=23/2-21/2 u 
J=55/2-53/2 a 
J=49-48 (003) [-3] 
J=49-48 (003) [1] 
J=15-14 


131267.422 
131285.894 
131286.219 
132179.132 
132381.516 


CH 2 CHCN 

c 5 h 2 n 1/2 

HC5N 

C 4 H 

U132378 


Jk+,K- = 14o,14-13o,13 

J=55/2-53/2 b 
J=49-48 (003) [+3] 
J=14-13 1/7 



Note. — Tabic A2 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 



data set has been compared with a complete model gen- 
erated from the basic physical conditions and chemical 
abundances deduced in the previous papers of this series 
(P04, P05, and PC06). Attempts have been made to 
search for new molecular species, not published in those 
works, but so far the only new finding is CH 2 CHCN in 
the outer Cold Circumstellar Shell of CRL 618, with a 
chemical formation route that must be similar to those 
of CH 3 CN and CH 3 CCH. The resulting morphological, 
physical and chemical picture of CRL 618 that has em- 
anated from this study is unprecendently precise thanks 
to the wide range of energy levels in the detected molec- 
ular lines that have subsequently been analyzed to probe 
the object. We have been able to provide details much 
beyond the angular resolution of the observations. Only 
intcrferometric observations could have provided a sim- 
ilar level of detail. However present-day interferometry 
could have concentrated on only a few lines and so a 
lot of information on the chemical composition would be 
missed. The present model can now be easily extended 



to the frequency range 280-360 GHz, for which another 
survey is almost completed with the Caltech Submillimc- 
ter Observatory. The Herschel Space Observatory will 
provide the opportunity to access higher frequencies into 
the submillimetcr range with an unprecedented sensitiv- 
ity and limited dilution compared with ISO. This will 
allow to access new molecular species and to complete 
the study of the continuum emission in this source. The 
advent of EVLA and ALMA will allow direct mapping of 
the structures delineated in this work and the possibility 
of a nearly continuous coverage from 300 MHz to 950 
GHz. The large collecting area should allow to extend 
our view of the chemical system. 

We acknowledge the support of the IRAM-30m staff 
during the long completion of the line survey. This 
work has also been supported by Spanish DGES and 
PNIE grants ESP2002-01627, AYA2002-10113-E and 
AYA2003-02785-E. 



APPENDIX 
ONLINE MATERIALS 

Tables Al, A2 and A3 contain the identified features in the CRL 618 millimeter wave survey in the following 
frequency ranges: 80.3-115.5 GHz, 131.2-179.1 GHz, and the 204.5-276.0 GHz. Also included are online only figures 
showing the millimeter spectra. 
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TABLE 3 

Identified features in the 204.5-276.0 GHz range in the CRL 618 survey presented in this paper. 



i/ (MHz) 


Molecule 


Transition/rJ-state 


i/ (MHz) 


Molecule 


Transition / w-state 


204788.938 


c-C 3 H 2 


Jk+,K- = 22,0-31,0 


204970.268 


HC 5 N 


J=77-76 


205003.000 


U205003 


205018.109 


CH 3 CCH 


J=12-ll K=4 


205045.497 


CH 3 CCH 


J=12-ll K=3 


205065.067 


CH 3 CCH 


J=12-ll K=2 


205076.813 


CHsCCH 


J=12-ll K=l 


205307.362 


HC 5 N 


J=77-76 (001) [-1] 


205308.360 


HC 5 N 


J=77-76 (010) [-1] 


205385.073 


HC 5 N 


J=77-76 (010) [+1] 



Note. — Table A3 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 
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TABLE 4 

CYANOPOLYYNE SPECIES INCLUDED IN THE MODEL AND J up AND/OR. E U p RANGE IN WHICH THEY APPEAR IN THE DATA. We SEPARATE 
ISOTOPIC SUBSTITUTIONS AND VIBRATIONALLY EXCITED STATES (THEIR VIBRATIONAL ENERGY IS ALSO GIVEN) OF THOSE SPECIES FOR WHICH 

SEVERAL OF THEM HAVE BEEN DETECTED. 



Molecule and 




1 

«J up 


Kangc trot 


lotal iNum. 


rvi /rup mi 
|AJ/|HG3l\IJ 


rvi /rim mi 


rvi /run mi 


vib. state 


(cm x ) 


range 


(cm J 


of labels 


in SEE 


in HVW 


in CCS 


HC 3 N (0000) 




9-30 


15-205 


17 


0.4092 


0.2325 


0.9905 


HC3N (0001) 


223 


9-30 


15-205 


34 


0.2412 


0.1859 


0.0095 


HC3N (0002) 


446 


9-30 


15-205 


51 


0.1065 


0.1115 




HC3N (0010) 


499 


9-30 


15-205 


34 


0.0531 


0.0598 




HC3N (0100) 


663 


9-30 


15-205 


34 


0.0216 


0.0304 




HC3N (0003) 


669 


9-30 


15-205 


68 


0.0418 


0.0594 




HC3N (0011) 


721 


9-30 


15-205 


68 


0.0313 






HC3N (1000) 


880 


9-30 


15-205 


17 


0.0033 






HC3N (0101) 


886 


9-30 


15-205 


68 


0.0127 






HC3N (0004) 


892 


9-30 


15-205 


77 


0.0154 






HC3N (0012) 


944 


9-30 


15-205 


100 


0.0138 






HC3N (0020) 


998 


9-30 


15-205 


49 


0.0052 






HC3N (1001) 


1103 


11-30 


24-205 


27 


0.0019 






HC3N (0102) 


1109 


11-30 


24-205 


88 


0.0056 






HC3N (0005) 


1115 


11-30 


24-205 


66 


0.0054 






HC3N (0110) 


1162 


11-30 


24-205 


62 


0.0029 






HC3N (0013) 


1168 


12-30 


29-205 


98 


0.0055 






HC3N (0021) 


1221 


15-30 


46-205 


64 


0.0031 






HC3N (1002) 


1326 


15-30 


46-205 


26 


0.0009 






HC 3 N (0103) 


1332 


15-30 


46-205 


90 


0.0023 






HC3N (0006) 


1338 


15-30 


46-205 


22 


0.0019 






HC 3 N (1010) 


1379 


16-30 


53-205 


9 


0.0004 




_ 


HC3N (0111) 


1385 


17-30 


60-205 


57 


0.0017 






HC3N (0014) 


1431 


17-30 


60-205 


56 


0.0016 


- 


- 


HC 3 N (0022) 


1464 


19-30 


75-205 


55 


0.0012 


- 


- 


HC 3 N (0030) 


1497 


23-30 


111-205 


20 


0.0005 


- 


- 


HC3N (1100) 


1543 


23-27 


111-165 


5 


0.0002 


- 


- 


HC3N (0201) 


1549 


23-30 


111-205 


19 


0.0005 


- 


- 


HC3N (1003) 


1549 


23-29 


111-191 


14 


0.0003 


- 


- 


rlU3JN (U1U4) 


1555 


23-29 


111-191 


21 


0.0008 






Up AT /I f!1 1 ^ 

riL/3lN ^iUii ) 


1 fim 
lOUz 


z4-oU 


1 01 one; 


Q 

o 


n nnnn 
U.UUUo 






HPoN (Tl1 1 9\ 


1608 


26-30 


14? 905 


21 


0.0007 






hc-jN ?nm 5^ 


1608 


26 


142-154 


1 


0.0007 






H 13 CCCN fOOOO"! 

i-J- V_/ V_ J V_/± ^ \\J\J\J\J ) 




10-31 


20-218 


18 


0.0273 


0.0155 


0.0248 


HC 13 CCN (0000) 




9-30 


15-205 


17 


0.0273 


0.0155 


0.0248 


HCC 13 C!N fOOOO^ 




9-30 


15-205 


17 


0.0273 


0.0155 


0.0248 


tj 13 r^rr^w /nnm \ 
ri OOOIN (UUUi ) 


ZZo 


in q 1 


on 01 q 
zU-Zlo 


OO 


n ni fii 






HC CCN (0001) 


223 


9-30 


15-205 


34 


0.0161 


- 


- 


HCC 13 CN (0001) 


223 


9-30 


15-205 


35 


0.0161 


- 


- 


H 13 CCCN (0002) 


446 


16-30 


53-205 


8 


0.0071 


- 


- 


HC 13 CCN (0002) 


446 


17-26 


60-154 


9 


0.0071 


- 


- 


HCC 13 CN (0002) 


446 


16-28 


53-178 


13 


0.0071 


- 


- 


H 13 CCCN (0010) 


499 


12-30 


29-205 


26 


0.0035 


- 


- 


HC 13 CCN (0010) 


499 


10-30 


20-205 


32 


0.0035 






HCC 13 CN (0010) 


499 


10-30 


20-205 


29 


0.0035 






H 13 CCCN (0003) 


669 






1 


0.0028 






HCCNC 




10-29 




15 


20-191 




0.0250 


TOTAL HC 3 N 








1736 








HC5N v=0 




31-104 


39-485 


55 


0.0132 




0.2794 


HC 5 N u V L 


105 


31-104 


39-485 


118 


0.0149 




0.0451 


HC 5 N 21/11 


210 


31-104 


39-485 


147 


0.0126 




0.0055 


HC 5 N 1/10 


230 


31-104 


39-485 


105 


0.0075 




0.0023 


HC 5 N 3i/ii 


315 


31-95 


39-402 


152 


0.0094 




0.0006 


HC 7 N u=0 




72-102 


102-198 


26 






0.0503 


HC7N I/15 


62 


79 


113-116 


1 






0.0052 


TOTAL HC 2 „+iN 








2340 
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TABLE 5 

Other molecules included in the model, presented as in Table [4] 



Molecule and 






Range E ro t 


Total Num. 


[Xl/fHC^Nl 

L J / L " J 


[Xl/fHC-iNl 


[Xl/fHC^Nl 


vib. state 


v / 


range 


(cm- 1 ) 


of labels 


in SEE 


in HVW 


in CCS 


CO 




1-2 


0-12 


2 


- 


667 


100 


13 co 




1-2 


0-11 


2 


2.5 


- 


2.5 


c 17 o 




1-2 


0-11 


2 


0.3 


- 


0.3 


c 18 o 




1-2 


0-11 


2 


0.2 


_ 


0.2 


HCO+ 




1-3 


0-18 


3 


0.1 


0.2 


0.1 


H 13 CO+ 




1-3 


0-18 


3 






0.003 


HOC+ 








1 








HON 




1-3 


0-18 


3 


2.0 


2.0 


1.5 


HON v2 


712 


1-3 


0-18 


5 








H 13 CN 




1-3 


3-17 


3 


0.133 


0.133 


0.033 


H 13 CN vl 

LA. 1 > 1/ _ 




1-3 




5 








HNC 




1-3 


0-18 


_• 


0.2 


0.2 


0.167 


fTNP 1/9 

111 N V_' IS 


477 


1-3 


3-18 


2 








TTT\Tl3p 




1-0 


n 1 s 





U.UIO 


n 01 


n nru 


ri2 w 






n 64 


n 






n 1 95 

O. lij 


SiO 




2-6 


1-30 


4 






0.013 


cs 




2-5 


2-35 


3 






0.125 


XT™ tr+ 

IN 2^1 




1 1 
1- 1 


U-o 


1 
1 






n 09 £ 


C-U3H.2 






9 1^/1 


fi7 






U.oD 


IvlglN 




7 1/1 


O-OO 


1 1\ 






n in 


0x130 IN 




0-10 


q qnn 
y-ouu 


OO 


U.UD / 


n 1 nn 

U. 1UU 


u.uoo 


V_y A ± J V_y ± J_ 




5-16 


6-315 


65 


1.0 


2.0 


0.200 


CH2CHCN 








120 






0.067 


CN 




1-2 


0-11 


26 






4.0 






1 9 










n 9p» 
u.zo 


CCH 




1-3 


0-18 


17 






0.5 






9-22 


1 4 77 


96 






nils 


C-C3H 






0-41 


25 






0.050 


1-C3H 






5-52 


58 






0.017 


C 4 H 




9-29 


14-142 


30 






1.0 


C 4 H u 7 


134 


9-22 


12-86 


31 








C 5 H 






23-100 


38 






0.005 


H2C3 








2 






<0.005 


TOTAL no-HC 2n +iN 








666 








H recombination 








44 








He recombination 








3 








U-lines 








74 








TOTAL LABELS 








3127 
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TABLE 6 

Unidentified features in the 80-276 GHz IRAM-30M survey of CRL 618. Frequency given in MHz. Area (A) given in K-kms" 1 

AND HPLW (W) GIVEN IN KMS" 1 . THE a OF THE FEATURE HAS BEEN ESTABLISHED FROM THE FOLLOWING AVERAGE RMS OF THE DATA: 4 
MK AT 3 MM. 8 mK AT 2 MM. 11 mK FROM 204 TO 240 GHz, AND 14 MK ABOVE 240 GHz. Only 74 UNIDENTIFIED FEATURES ARISE ABOVE 

THE 3<T LIMIT. 



Freq. 


A (W) cr 


Freq. 


A (W) a 


MHz 


K km s^ 1 


MHz 


K km s- 1 




(km s" 1 ) 




(km s" 1 ) 



80559 


0.129 (18) 3a 


80618 


0.069 (7) 3a 


83659 


0.380 (24) 4cr 


85136 


0.295 (16) 5a 


87013 


0.160 (17) 3a 


89120 


0.124 (9) 4a 


89600 1 


0.239 (27) 3a 


89820 


0.124 (9) 4a 


91156 


-0.320 (20) 4a 


91238 2 


-0.082 (5) 4a 


91383 


0.389 (21) 5a 


91441 2 


0.100 (9) 3a 


92414 2 


-0 105 (4) 7a 


92814 3 


-0 147 (7) 5a 


95667 


0.289 (23) 3a 


97637 


0.447 (48) 3a 


97728 


0.167 (17) 3a 


97772 


0.246 (22) 3a 


98822 


0.301 (27) 3cr 


98856 


0.392 (34) 3a 


100265 


0.304 (18) 4a 


100540 


0.530 (37) 4a 


101372 


-0.115 (5) 6a 


101751 


0.212 (14) 4a 


102222 


0.261 (18) 4a 


103426 


0.160 (17) 3a 


103910 4 


0.237 (19) 3a 


107121 5 


0.281 (24) 3a 


107275 


0.310 (25) 3a 


107620 


0.377 (31) 3a 


115426 


0.675 (19) 9a 


132378 


0.22/-0.49 (9/27) 4a 


132679 


0.735 (52) 3a 


134239 


0.248 (11) 3a 


134525 


0.488 (28) 3a 


137291 


-0.121 (4) 4a 


137720 


0.446 (15) 4a 


141811 


0.117 (8) 3a 


142293 


0.165 (9) 3a 


146210 


0.445 (17) 3a 


150349 


-0.054/-0.055 (3/6) 3a 


150620 


0.470 (42) 3ct 


160670 


0.835 (27) 4(t 


160952 6 


3.669 (56) 9a 


170642 


0.847 (35) 3a 
-0.389 (15) 3a 


172998 


-0.092 (4) 3a 


173106 


173572 


-0.130 (6) 3a 


175840 


0.057 (2) 3a 


205003 7 




218042 


1.144 (17) 6a 


218207 


1.640 (17) 9a 


219270 


0.575 (18) 3a 


220431 


0.556 (17) 3a 


221366 


2.132 (59) 3a 


221980 7 


3.612 (28) 11a 


228507 


0.105 (3) 3a 


234936 


0.693 (11) 5a 


235597 


0.422 (16) 3a 


235606 


0.152 (6) 3a 


241337 8 


3.799 (45) 7a 


241388 s 


1.175 (19) 5a 


241498 


0.913 (15) 5a 


242642 


0.178 (5) 3a 


245120 


0.690 (17) 3a 


248647 


0.139 (4) 3a 


249111 


0.487 (17) 3a 


256428 


0.238 (4) 5a 


258314 


0.811 (15) 4a 


260295 


0.293 (6) 4a 


260426 


0.233 (5) 3a 


263475 


0.328 (8) 3a 


263840 


0.186 (4) 3a 


272988 


1.560 (11) 10a 



1 (id) HC la CCCCN J=34-33 x (cm) only one channel 

3 (id) C 3 H 2 4 (id) H 56 7 

5 (id) HC 5 N 3i/io(l+) J=40-39 6 (id) CH 2 CN J=8-7 

7 (bl) CH3CCH 8 (id) CH 2 CN J=12-ll 



id: possible identification; cm: comment; bl: existing blending 
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Fig. 1. — Results of the CRL 618 IRAM-30m survey in the frequency range 80.0-116.5 GHz (3 mm window) compared to the model 
devepoloped according to the results found in the series of papers P04, P05 and PC06. We include nesting zooms on intervals covering 4 
GHz and 0.5 GHz respectively. The data are shown in terms of with the continuum level substracted. The labels next to the arrows 
in the upper panel indicate the J up number of the corresponding HC3N or HC5N transition. The complete set of detaliled figures at 0.5 
GHz intervals with line identifications (similar to the lower panel), is available as on-line material only. 
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Fig. 2. — Same as previous figure but for the frequency range 130.5-179.0 GHz (2 mm window). The first zoom shows 8 GHz instead of 
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Fig. 3. — Same as previous figure but for the frequency range 204.5-245.0 GHz (lower part of the 1.3 mm window). 
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Fig. 4. — Same as previous figure but for the frequency range 245.0-275.5 GHz (upper part of the 1.3 mm window). 



